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Abstract 

We study gravitational waves propagating through an anisotropic 
Bianchi I dust-filled universe (containing the Einstein-de-Sitter uni- 
verse as a special case). The waves are modeled as small perturbations 
of this background cosmological model and we choose a family of null 
hypersurfaces in this space-time to act as the histories of the wave- 
fronts of the radiation. We find that the perturbations we generate 
can describe pure gravitational radiation if and only if the null hy- 
persurfaces are shear-free. We calculate the gauge-invariant small 
perturbations explicitly in this case. How these differ from the corre- 
sponding perturbations when the background space-time is isotropic 
is clearly exhibited. 
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1 Introduction 



In a recent paper shear-free gravitational radiation analogous to the 
Bateman waves of electromagnetic theory was studied propagating through 
isotropic cosmological models. The radiation was represented as a pertur- 
bation of the cosmological model. It was analysed using the Ellis-Bruni 
gauge-invariant and covariant perturbation formalism. Explicit solutions 
of the linear differential equations for the gauge-invariant small quantities 
(the Ellis-Bruni variables) were derived. These variables divide naturally 
into so-called vector quantities and tensor quantities. Examples of the vec- 
tor quantities are the spatial gradients of the proper density of matter, of 
the isotropic pressure, of the expansion scalar of the matter world-lines, as 
well as the vorticity vector associated with the matter world-lines and the 
heat flow vector of the matter distribution. The tensor quantities include 
the shear of the matter world-lines and the anisotropic stress of the matter 
distribution. For perturbations describing gravitational waves the important 
variables needed are the tensor variables []J, ||. If the cosmological model is 
isotropic then the shear of the matter world-lines and the anisotropic stress 
of the matter distribution are Ellis-Bruni variables. However if the cosmolog- 
ical model is anisotropic then the matter world-lines have shear and so the 
shear tensor is not (in its entirety at least) available to us as an Ellis-Bruni 
variable. This has important implications therefore for the study of gravita- 
tional wave propagation in an anisotropic universe. In this paper we show 
how to circumvent this difficulty In addition we discover what influence the 
anisotropy in the universe has on the gravitational waves. 

We are concerned here with deriving perturbations of an anisotropic cos- 
mological model which can be interpreted as representing shear-free gravita- 
tional radiation. For definiteness we choose a Bianchi I dust-filled universe 
(containing the Einstein-de Sitter universe as a special case) and a simple 
family of null hypersurfaces in this space-time to act as the histories of the 
wave fronts of the radiation. To begin with these null hypersurfaces can have 
shear but we show that the perturbations we generate can only represent 
gravitational waves (with the gauge-invariant part of the perturbed Weyl 
tensor having the algebraic structure one associates with pure gravitational 
radiation) if the null hypersurfaces are shear-free. This involves a special- 
isation of the anisotropic universe. The perturbed quantities are explicitly 
calculated in this case and all are derived from a single complex-valued func- 
tion whose dependence on one complex variable is analytic and arbitrary (a 
feature always associated with shear-free gravitational radiation) and whose 
dependence on two real variables is determined by a wave equation which 
can easily be integrated. If the background cosmological model is isotropic 



2 



(Einstein-de Sitter) then the wave equation is second order whereas if the 
background is anisotropic the wave equation is first order. This difference is 
directly attributable to the presence of shear in the matter world-lines in the 
anisotropic background. 

Density perturbations of anisotropic Bianchi I universes have been thor- 
oughly examined in the Ellis-Bruni formalism by Dunsby [|J. Our study of 
gravitational wave perturbations of a Bianchi I universe is complementary to 
this work and has in part been inspired by it. 

The outline of the paper is as follows: The notation used throughout 
this paper and some important basic equations are given in Appendix A 
and referred to where appropriate. In Sec. II we describe the Bianchi I 
dust-filled universe which will play the role of a "background" space-time 
for our small perturbations. In Sec. Ill we identify the Ellis-Bruni gauge- 
invariant variables we shall use and give the equations satisfied by these 
variables. Then, in Sec. IV we require these Ellis-Bruni variables to have an 
arbitrary dependence on a function (in order to have the possibility of having 
gravitational waves with an arbitrary profile) and insert this dependence into 
the equations. The calculation of the perturbed quantities when the null 
hypersurfaces, described in the second paragraph of this introduction, are 
shear-free is also given in Sec. IV and the analogous calculations for the 
shearing hypersurfaces (also referred to above) are given in Appendix B. 
Finally, Sec. V is a discussion in which we summarise and comment on the 
results obtained in Sec. IV and in Appendix B. 

2 An Anisotropic Space— Time 

We are interested in constructing models of gravitational waves propagat- 
ing through an anisotropic universe. For clarity we choose a specific Bianchi 
type I dust -filled universe which is the dust -filled counterpart of the Kasner 
universe (see, for example, 0) . The space-time model has line-element 

ds 2 = A 2 (t) dx 2 + B 2 (t) dy 2 + C 2 (t) dz 2 - dt 2 , (2.1) 

with A(t) = (t-k)P(t + k)-3~ p ,B(t) = (t-k)«(t + k)*-*,C(t) = (t-k) r (t + 
k)s~ r , where p, q, r are constants satisfying p + q + r = 1 and p 2 + q 2 + r 2 = 1 
and k is a constant (if k = the line-element ( |2.1|) becomes the Einstein- 
de-Sitter line-element). 

The world-lines of the dust-particles are the time-like geodesic integral 
curves of the vector field u a d/dx a = d/dt (thus u a = 51 since we shall label 
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the coordinates x 1 = x, x 2 = y, x 3 = z and x 4 = t). The proper density is 

AB AC BC _ 4 
fi ~AB + AC + BC~ 3 (t 2 - W) ' (2 ' 2) 

where a dot indicates differentiation with respect to t. For future purposes 
it is useful to note that the field equations satisfied by the functions A, B, C 
given above are 

), (2.3) 
), (2.4) 



BC 


B 


C 


BC 


+ B 


+ c 


AC 
AC 


A 
+ A 


c 
+ c 


AB 
AB 


A 
+ A 


B 
+ B 



3. (2.5) 

The geodesic world-lines of the dust particles are twist-free (vorticity tensor 
^ab = 0) and have expansion scalar 

„ A B C 2t , . 

e -A + B + c = ¥^¥' (2 - 6) 

and shear tensor 

O-ab = 0"(l)n(i) a n ( i) 6 + <T(2)Tl(2) a n(2)6 + 0-(3)"(3)a"(3)6 , (2-7) 

where n (1)a = A 5^ , n (2 ) a = B5 2 al n {3)a = C 5% and 



with A, 5, C given above and given by ( |2.(j| ). The vorticity and shear 
tensors and the expansion scalar are defined in Appendix A. 

We note that a a b = aba satisfies a ab = and, in addition, a a a = (this 
is equivalent to am + a^) + o"(3) = which follows for example, from (|2.6|) 
and (|2.8|) ). Following Ellis it is convenient to define a projected covariant 
derivative operating on tensors which are orthogonal to u a on all of their 
indices. In particular, if for example T a b c is such a tensor, then its projected 
covariant derivative, denoted by ^V^T^ , is defined as 

{Z) V d T abc = h p d hlhlh s c T qrs , p , (2.9) 

where h% = 5% + u a Uf, is the projection tensor and the semi-colon indicates 
covariant differentiation. When this derivative operates on a scalar function 
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/ it naturally specialises to the partial derivative projected orthogonal to u a . 
Thus 

{3) ^af = h b a f, b , (2.10) 
with a comma denoting the partial derivative. We note that 

(3) V a /i bc = 0, (2.11) 

and thus the projected covariant derivative commutes with the operations of 
lowering and raising indices on tensors orthogonal to u a using h ab and h ab 
respectively. 

In the exact Bianchi type I universe the projected covariant derivative of 
any scalar function vanishes. Thus we have 

(3) Va # = 0, ( 3 )V ay u = 0. (2.12) 

Also 



(3) V fe a (a) =0, (2.13) 



for a = 1, 2, 3, and hence 



( 3 )V 6 a 2 = 0, (2.14) 

where a 2 = \ a ab & ab = ofn + of 2 ) + cr%) ■ ^ follows from ( pT7[ ) , ( [2.1 3D and the 
fact that {3) V c n a {a) = with a = 1, 2, 3 that 

( 3 )V c a ab = 0. (2.15) 

With 9 and a a b given by (p.6|)and (|2.7| ) respectively, it is easily shown using 
0, (PD, and (p]) that 

a ab + 8a ab = 0. (2.16) 

Here and throughout a dot indicates covariant differentiation in the direction 
of u a . Thus & ab = a ab . c u a . 

With respect to the 4-velocity field defined after ( [2.1| ) the Weyl tensor 
with components C abc d, is decomposed into an "electric part" and a "magnetic 
part" given by 

E ab = C apbq up , H ab = *C apbq u p u q . (2.17) 

Here *C apbq is the dual of the Weyl tensor and it is defined by *C apbq = 
lVa P rs C rsbq , where 7] abcd = yj=g ' t abcd with g = det (g ab ) and e abcd is the 
Levi-Civita permutation symbol. The electric part of the Weyl tensor of the 
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anisotropic universe with line element ( |2.1| ) is non-zero and can be written 
in terms of the shear a ab of the matter world-lines as 

2 1 
E a b = g cr 2 hab ~ cr ac a c b + -9 a ab . (2.18) 



This follows from the shear-propagation equation (see equation ( |3.5[ ) below) 
in this space—time using the facts that the matter world-lines in this case 
are geodesic, twist-free with shear satisfying ( |2.16| ) and the dust matter 
distribution is stress-free. The magnetic part of the Weyl tensor vanishes in 
this universe and so we have 

H ab = 0. (2.19) 

In the rest of this paper we wish to construct perturbations of this 
anisotropic universe which we can interpret as gravitational waves propa- 
gating through the universe. The histories of the wave-fronts will be null 
hypersurfaces in this anisotropic (background) space-time. To be specific we 
shall choose the null hypersurfaces to have equations 

<j) (x, t) = x — T(t) = constant , (2.20) 

with dT/dt = A" 1 . It is easy to see that these hypersurfaces are null with 
respect to the metric g a b given via the line-element ([2.1|). Thus 



g ab <f>, a <t>,b = 0- (2-21) 

The integral curves of the vector field g ab 4>,b are the null geodesic gener- 
ators of the hypersurfaces ( |2.20| ). In terms of a convenient null tetrad in the 
anisotropic background given by the 1-forms, 

k a dx a = Adx — dt, l a dx a = —-(Adx + dt) , 

1 1 
m a dx a = -j={Bdy + iC dz) , m a dx a = -j= (Bdy - iCdz) , (2.22) 

the complex shear of the null geodesies tangent to 4>,a is 

a = <f ) , a ., b m a m b = ^ ( | - = -L (a {2) - a (3) ) , (2.23) 

where we have used (|2.8|) in the last equality. We see that the null hyper- 
surfaces are shear-free if and only if B = C . With the explicit expressions 
for A, B, C given following ( [2.1|) we see that this corresponds to putting the 
constant r = q and choosing q and p such that p + 2 q = 1 and p 2 + 2 q 2 = 1. 
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3 Equations for Perturbations 



We wish to study gravitational wave propagation in the anisotropic uni- 
verse described in Section II. We use the Ellis-Bruni approach 0] which 
involves working with small gauge-invariant and covariant quantities which 
vanish in the background rather than small perturbations of the background 
metric. The gauge-invariant variables used depend on the type of perturba- 
tions to be studied and also on the background space-time. For example, 
E ab defined by fl2.17[ ) is an Ellis-Bruni variable for gravitational wave propa- 
gation in isotropic cosmologies U , but in the present context its background 
value is non-zero (it is given by ( |2.18| )) and hence it is not an Ellis-Bruni 
variable for our purposes. However, it is possible to extract a gauge-invariant 
part E ab of E ab given by (see equations ( |3.13| ) and Q3.14D below) 



E a b ■= - n ab - (& ab + 6a ah ). (3.1) 

Here n ab is the small anisotropic stress in the perturbed matter distribution 
(see Appendix A). For our case the important Ellis-Bruni variables are E ab , 
H a b, Tr ab , a ab + 6a ab and ^V c <J a b- We set all other Ellis-Bruni variables, 
including q a , u ab , h b c 9 yb , h b c fi tb , h b c p^ b {p is the isotropic pressure) and u a , 
equal to zero since we found in JTJ that it is tensor, not vector, quantities 
that describe gravitational wave perturbations (the vector quantities describe 
other types of perturbations such as inhomogeneities in the matter density). 
We now outline the equations satisfied by these variables. 
When the Ricci identities, 

U a ;dc — U a]C d = R a bcdU b , (3.2) 

where R a bcd is the Riemann curvature tensor, are projected along u a and 
orthogonal to u a (using the projection tensor) we obtain Raychaudhuri's 
equation, 

+ \e 2 - u a , a + 2 (a 2 - lu 2 ) + \{n + 3p) = , (3.3) 
(here u 2 = \ uj ab ui ab ) the vorticity propagation equation, 

h a b u b +^6u a = a\ u b + l - V abcd u b u c , d , (3.4) 

(where cu b := | rf bcd u b u cd is the vorticity vector) the shear propagation equa- 
tion 



2 

K K (cr fg - U(f ;g )) -ii a u b + u a uj b + cr af a f b + - 9 a ab 

+h ab {-- u 2 --a 2 + - u c -A - - 7c ab + E ab = , (3.5) 
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the (0,u)- field equation (this terminology is explained in |J), 

~ hi 9> b - h a b a bc , d h d c - r?^ u c (cu d .j + 2uj d u f ) = q\ (3.6) 

the divergence of vorticity equation, 

u a , b h b a = u a u a , (3.7) 

and the magnetic part of the Weyl tensor, 

H ah = 2 u {a u b) - ti a h s b {u {t 9 ' c + a {t 9 ' c ) rj s)fgc u f . (3.8) 

Next we consider the conservation equation, 

T ab . b = 0, (3.9) 

where T ab is the energy- momentum-stress tensor (see Appendix A). This 
equation projected orthogonal to u a and along u a respectively gives the equa- 
tions of motion of matter, 

(jm + p)u a + h ac (p )C + 7T b c , b + q c ) + [u ab + a ab + U h ab ^j q b = 0, (3. 10) 

and the energy conservation equation, 

fi + e(fx + P ) + 7f ab a ab + q a . a + u a q a = . (3.11) 



Specialising equations (|3.3|) - (|3.8|) and ( |3.10|) -( j3.11| ) to the problem at hand 



with perturbed quantities listed after ( |3.1| ) and q a , u ab , h b c 9 b , h b c fi !b , h b p <b 

and u a all put equal to zero, we obtain the following equations: From Ray- 
chaudhuri's equation, 

9 + \e 2 + 2a 2 + ]-(fi + 3p) = . (3.12) 
From the shear propagation equation, 

E a b = ~ TT a b ~ 0~ ab ~ &af ^ b - - 9 O ab + - ff 2 h ab . (3.13) 

With E ab given by ( |3.1[ ) we can write this equation as 

Eab = E ab - a af a f b + ^ 9 a ab + ^ a 2 h ab . (3.14) 



We now see explicitly from equations ( |2.19|) and (|3.14j) that E ab = in the 



anisotropic background. This justifies the definition of E ab in equation ( |3.1[ ). 



From the (0,v)- field equation 



(3) Va a a6 = 0. (3.15) 

From the magnetic part of the Weyl tensor 

H ab = -hl a ^\/ c aP b)Vsfpc u f . (3.16) 

From the equations of motion of matter 

h ac 7r b c , b = 0. (3.17) 

Finally the energy conservation equation gives 

fi + 9(fi + p) + n ab a ab = . (3.18) 

The remaining two equations derived from the Ricci identities (the vorticity 
propagation equation and the divergence of vorticity equation) are identically 
satisfied. We note that equations (|3.12|) and ( |3.18|) are not expressed in terms 
of the Ellis-Bruni variables we wish to work with and so are not immediately 
useful. However, taking the projected covariant derivative of these equations 
and retaining only linear terms yields: the projected covariant derivative of 
Raychaudhuri 's equation, 

or ab {3) V c cr ab = , (3.19) 
and the projected covariant derivative of the energy conservation equation, 

a ab ^V c ir ab = . (3.20) 
When the Bianchi identities, written conveniently in the form 

C abcd d = R c[a;b] _ g c[a R ;b] ^ ^1) 

where R ca := R cba b are the components of the Ricci tensor and R := R° c 
is the Ricci scalar, are projected along u a and orthogonal to u a we obtain 
equations for E ab and H ab which are analogous to Maxwell's equations. The 
general forms of these equations are lengthy and to keep this section to a 
reasonable length we give them in Appendix A. In terms of the perturbed 
quantities given following ( |3.1| ) and with q a , u ab , h b 9 b , h b fi )b , h b c p b and u a 
all vanishing these equations are: the div-E equation (using ( |3.17| )), 

( 3 ) V a E ab = a a f ®V a a b f + V ba ™ u a a\ H qd , (3.22) 
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the div-H equation, 

(3) Va H ab = - V ba ™ u a a d p (E qd + 1 vr^ , (3.23) 

the H -equation, 

H bt - h {b n t)rsd u r E a s . d - r] rsd{t u r (3) V d a b)s a fs - rf sd( * u r a b)f (3) V d a fs 
~0 ( (3) V d a (6 s ) rf )rsd u r - 3 H^s a b)s + /i 6 * #* a dp + 

= - ~ V [b ra d vr t)a;d « r + i r/ 6 ra(i u *) u r o~ cd 7l\ , (3.24) 



and the E -equation, 

-V - la + -fij a +6cr f a J s a - h a f cr J a da 

= -Ebt _ 2 ^ + ^ ^ _ d E bt + {& b f + 6 a b f ) a ft 
3 

+3 #* c a b ) c - /»" E* a dp + \ h bt ** a dp - U bt - \a< b A 

--9n bt . (3.25) 
6 

We note that the right-hand side of equation ( |3.25| ) is expressed in terms 
of Ellis-Bruni variables and thus the left hand side must be an Ellis-Bruni 
variable which we denote by 

W bt := Q 6 - 2 a 2 + - /x) a bt + 3 a b > ^ a 5 * - fc w <x /d <x ad . (3.26) 

Calculating the projected covariant derivative of W bt gives 

&V c W bt = -3a 2 ^V c a bt + 3^V c (a b f a f s a st ) 

-h bt ®V c (a a f o-f d a da ) , (3.27) 

where we have used the background values of 9, \i and a ab given in the 
previous section to simplify the coefficient of ^\/ c a bt . Using ( |3.25| ) and 
(|3.1| ) we can express W tb in terms of E bt , H bt and n bt as 



~ bt 2 1 1 1 

W bt = —E - - 6E bt - -7r bt - - 6ir bt - - h bt ir dp a dv 
3 2 3 6 p 

± h bt E dp a dp + l - aOf n t)f + E<-* f a b)f - h<£ ^ rsd u r H a s , d . (3.28) 
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We shall use this expression later in the left-hand side of ( |3.27| ). 

Next we examine Eq. (|3.23| ). Substituting for H ab from (|3.16| ) gives 

V frsd u r {-\ pV d (7* f ) !/ + tt» ( r /8 Wv d a B .}-^- ffre (»V ( ,<7* l 

( {3) V d a a s ) a -a a s ^E da + ^7r ad ^ =0 . (3.29) 



-if rsd u. 



Multiplying across by r\umq an d using 



u b 



h 

d<? s 



--^ (3) V d ^ - a bl ®V d a b 



(3.30) 



(since ^*V d a b s is orthogonal to u b ) and 



f(3)y 
I V r 



ab (3) Vq 



(3) 



m & q I j 



(3.31) 



we find that the div — H equation implies that 



(3) V b ( (3) V [m a fc g] 



+ ^{m [ E q ] a + - 7T g ] a 







(3.32) 



with the square brackets as always denoting skew-symmetrisation. The con- 
verse is also true. This is shown by multiplying (|3.32|) by rf smq and using 



(3)y a b _ (3)y J> — n v 

v m u q v q u m i/qpfm •* 



(3.33) 



which is easily derived from (|3.16|) . Thus we can replace ( |3.23|) by (|3.32|) 
without any loss of information. 

The projected covariant derivative and the covariant derivative in the 
direction of u a do not commute. The commutation relation satisfied by these 
derivatives when acting on <j ab follows from the Ricci identities and reads: 

(3) V c a ab - (( 3 ) V c a ab ) ■ = a\ ®V b a cs - a\ ^V s a cb + a\ ®V a a cs 

-a s b ^V s a ca + a s c {3) V s( T ab + le^V c a ab . (3.34) 



In the present context this important equation can be written solely in terms 
of gauge-invariant variables. Thus we shall use it in the form, 



(3) 



V c ( & ab + 9 a ab ) = ( (3) V C a ab ) ' - a f b (3) V/ a ac - a f a (3) V/ a bc 



+a f b {3) V a Of C + a f a (3) V b a fc + a f c (3) V/ a ab + - 9 (3) V c a ab . (3.35) 
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Using this equation and ( |3.1| ) we obtain the following equation for the pro- 
jected covariant derivative of E ab , 



(3) V c a ab ) ■ + a f b {3) V f cr ac + a f a (3) V/ a bc 



{3) V c E ab = ^V c 7r ab 
-o\ ®V a a fc - ® V 6 a fc - of c ^ V f a ab - ^ 9 (3) V c a ab . (3.36) 

Contracting this equation over a and c and making use of (|3.15|) and (|3.17|) 
gives 

W\7 a E ab = Q. (3.37) 

With H ab given by ( |3.16 ) it is easily shown that ( 3.22| ) is consistent with this 
equation. It also follows from ( |3.36| ), ( |3.19| ) and ( ^.20[ ) that 



a 



^V r E ab = 0. 



(3.38) 



The equations we will use are Q, (|3~Tj ), (|3~l5| ), fl5IE|) , ( gl9|) , 
(|33q) , ( |3gg ), |T27|), (ggg ), (|3~^ ), (|3~35| ), (|3TB1 ) and (|338| ). These equa- 
tions are not independent of each other. For example (|3.24j) is automatically 
satisfied and this can be seen as follows: First using ( |3.13| ) and ( |3.35| ) we 
show that 



rf rsd u r ^V d a b s ) ■ - 9 V trsd u r ^V d a b 



h aV trsd u r ^E a s --v\ 
- V trsd u r af s ( < 3 > V 6 a fd - W V, a b d + ^V d a b f - W V, a b d ) 



(3.39) 



Then substituting for ^V d cr b f — ^Vfcr b d from ( |3.33| ) in this equation we 
find that 



h^rf> 8d u r (E a s - l -Tt\ 



H bt -3a {t d H b)d + h bt a dp H dp -9H 



bt 



(3.40) 



If we replace E a s in this equation with the right-hand side of ( |3.14 ) the result 
is Eq. ( |3.24|) and thus the H — equation is automatically satisfied. We shall 
demonstrate the internal consistencies of the remaining equations by finding 
a solution which satisfies all the equations. 



4 Shear— Free Gravitational Waves 

We shall now look for solutions ^V c a ab , a ab + 9 a ab and 7r a6 of the equa- 
tions given in the previous section for which these variables depend upon 
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an arbitrary function. We expect that this dependence of perturbations will 
describe gravitational waves carrying arbitrary information. Specifically we 
assume that 

(3) V c a a6 = A cab F(<p) + G cab F'(<p) , 
& ab + 6a ab = B ab F{4>) + C ab F\4>) , (4.1) 

V"afe = n ab F(0) , 

where F is an arbitrary real-valued function of its argument <p{x a ) and F' 
denotes the derivative of F with respect to its argument. This idea of in- 
troducing arbitrary functions into solutions of Einsteins equations describing 
gravitational waves goes back to work by Trautman M. Its use in the con- 
text of gauge-invariant perturbations of cosmological models was initiated 
by Hogan and Ellis |3j and further developed by Hogan and O'Shea [T]]. We 
note that all of the quantities in ( |4.1| ) are orthogonal to u a and B ab , C ab and 
U a b are trace-free with respect to the background metric g ab (i.e. B a a = 0, 
C a a = and U a a = 0). 

When we substitute the first two equations in Eq. flj.l|) into the com- 
mutation relation ( |3.35| ) we, in effect, obtain integrability conditions to be 
satisfied by the right hand sides of these two equations. Remembering that 
F is an arbitrary function this substitution results in the following equations: 



where A c = h b c 4> tb and <p - 

B a b A c + ®V C C 6 = 

and 

(3) V c B ab = 

If we now let 

where s ab u a = and s a a 



C a b A c — G cab (j) , 



(4.2) 



- A cab (f) + G C ab — V b Gfac — <J a G fbc 

+V f b Gafc + V f a G bfc + (7 f c Gf ab + ^6 Gcab ,(4.3) 

Acab - CT f b Af a c - a f a Af bc + & f b A a f c 

+a f a A bf c + (T f cA f ab + -6Acab ■ (4.4) 

Cab = <pS a b , (4.5) 

= then it immediately follows from ( |4.2j ) that 

Gcab = Sab A c . (4.6) 
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As a result of ( |4.5| ) and (fO|) and also since 



(3)W 



and 



(4.7) 
(4.8) 



I now reads 
B a b A c — A cab < 



' + 4> u c + - 6 A c + o v c \ , 



(A c )' = (0,c)'+ 4>U C 



---^V c s ab -\6\ c s ab + s ab \ c -af h \ f s ac 

-<J f a \fS bc + (T f b \ a S fc + (T f a X b Sf c + -9\ c Sab -(4.9) 



Substituting Eqs. ( flip (with C ab and G cab replaced by ( |4.5| ) and Q4.6| ) re- 

spectively) into Eqs. tH), (EH) Q)> fUD' (H)> (UZD, fUD' 

( |8.32| ), ( |3.37| ) and ( |3.38| ) results in the following list of equations: From the 
(0, u) field equation, 

s ab X b = 0, (4.10) 
A b ab = . (4.11) 
From the magnetic part of the Weyl tensor, 

H a b = QabF + l ab F , 

where as always F' = dF /dip and 

q ab = - A cp {b r] a)fpc u f , 

Lb = ~ X c s p {b rj a)fpc u f . 
From the equations of motion of matter, 



nab „,a be _ 

■ h = U IT a bc 



rr 6 <j) >b = 

From the projected covariant derivative of Raychaudhuri 's equation, 

a ab A cab = 0, (4.17) 

a ab s ab = 0. (4.18) 
From the projected covariant derivative of the energy conservation equation, 

a ab ^V c Tl ab = , (4.19) 



(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 
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a ab U ab = 



From ( HOD 



From (13.3 



1 



£ ab = (^n ab -£ afe )F-0 Safc F\ 



(4.20) 
(4.21) 



qa J 



(4.22) 



and 



A m \ X b + X m (3) V fe s\ + A (3) V b X m - A b X b - A a (3) V b s b 



q I a q 

~S m ^ Xq (f) <J m Sq a -\- (j) <J g S ma . 



From fl3~37p 



(3) V aJ B afe = 0, 
B ab X a + 4> (3) V a s ab + s ab (3) V a = 
v 6 



s ab X b = 



From (ETggp 



From (gg?D 



(3) V,W 



6i 



a afc ( 3 )V c 5 ab = , 
a a6J B ab = . 

3 (_ ff 2 Ac te + (7 6 / ^/^ A/i + <J b f a st + ^^st A b f 



(4.23) 

(4.24) 
(4.25) 
(4.26) 

(4.27) 
(4.28) 



-h a n a 



+ 6{—a s + a f cr J s s + a fa s J s 



+a f s a st s b f-h bt a d a a a f s fd ) X C F' . 



(4.29) 



Finally, from Eq. ( ggg) using (P3)-(4.14), (P3) , (4.20), flOID and (4.28) 
we find that 



W 



bt 



■-ibt 2 w £ 6t 2 w (6 t }/ _ ( t b)f 
3 3 1 S 



{n 

(3) V(j ( A M 



A 



Ulb 



A bdt )}F + {(B bt -U bt + s bt 



+ ±0 s bt_ s ( t f a b) f) ^ + £ s bt _ A dbt ^ _ X d (3) Vd s bt 

3 

_ s i>t 0) v x d +- A tdb X d + -X l ®V d s db + - s db ^V d A* 



1 

2 

12 M 



V-A bdt X d 



I s ^(3) VdA 6 + l A 6(3) VdS dt }F / + 



{0 2 s w - A d A d s w }F" . 



(4.30) 
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Replacing W bt in ( |4.29| ) by the right-hand side of ( [4 . 3 Q| ) and equating the 
coefficients of F'" yields: 

(0 2 s bt - \ d \ d s bt ) X c = . (4.31) 

Assuming A c 7^ this implies that 

<P, a <P' a = 0. (4.32) 

Thus the hypersurfaces <fi (x a ) = constant in the background anisotropic 
cosmological model are null. The equations found here by equating the co- 
efficients of F, F' and F" are extremely complicated and we will work with 
them in a completely different way below. 

For the remainder of this section we seek to construct perturbations of 
the anisotropic background cosmological model which describe shear-free 
gravitational waves having the null hypersurfaces <fi(x a ) = constant as the 
histories of their wave fronts in the background space-time. For the function 
4>(x a ) in (|4.32| ) we shall take (p(x, t) given by fl2.20|) and in order to have 
these null hypersurfaces shear-free we see from ( |2.23| ) that we must have 
°"(2) = ct(3)- This is a relationship between the (non-vanishing) principal 
shears of the background matter distribution. One immediate consequence 
of this is that E a b is no longer the full gauge invariant part of E a b. To see 
this we note that in a general space-time with metric g a b and a preferred 
congruence of world-lines tangent to u a we can write the shear tensor a a b in 
terms of the principal shears <7(i) , at 2 ) , C(3) as 

&ab = C(i) 71(1) a 7i(i) b + CT(2) ™(2) a ™(2) 6 + C(3) "(3) a "-(3) b , (4.33) 

where <7(i) + o"( 2 ) + o"( 3 ) = and a , nm a and n^) a are the unit orthogonal 
(space-like) eigenvectors of a ab . Then letting m a = {n%% + infy)/^ and 



using ( |4.33| ) we can write 

(Tab ~ (T a f ° f b ~ g a 3 Kb^j m a m b = o ab m" m b , (4.34) 



1 ^ , 2 r> , \ „ f, / I 

3 

where ^ 

a ab m a m b = - (a {2 ) - cr {3) ) , (4.35) 

is now a gauge-invariant variable which vanishes in the background space- 
time in which am = cp). When <j a b is the perturbed shear m a is given its 
background value which, for our case, is given implicitly in ( |2.22| ). Thus we 
can give + | 9 its background value when (|4.34|) is calculated in first ap- 
proximation. It follows from this and ( |3.14| ) that now in first approximation 
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the full gauge invariant part of E a b is 

A 

J0 r , 



EL 



ab 



E a b + ~7 &rs ^ frf ™>a ™,b + ~T 0"rs ^ ™" ™a ?™& • 



.4 



(4.36) 



We shall consider pure pure gravity wave perturbations of the anisotropic 
background. We do this by requiring that the gauge-invariant perturbed 
"magnetic" and "electric" parts (H a b and E' ab respectively) of the Weyl tensor 
be type N in the Petrov classification with <j) ,a as degenerate principal null 
direction. We achieve this provided E' a 



ab ' 



= H, 



in light of ( 4.1C ) and ( 4.16 ), we should require 



B 



ab 







ab 







ab ' 



q ab 4>, b 



This means that, 



(4.37) 



with l ab , q ab given above in ( |4.13| ) and ( f4.14| ). The first of these allows us to 
write Eq. ( f4.25|) as 



(3)y 



-,ab\ 



which we can simplify to 



•s ab \b 







(4.38) 



(4.39) 



We note that it follows from the definition of l a b given in ( f4.14|) that the 
second of Q4.37j ) is identically satisfied. 

Since s ab , U ab and B ab are orthogonal to u a and <p' a and trace-free with 
respect to the metric tensor given via the line-element (|2.1| ) each have only 
two independent components, — s 22 = s 33 = a (x, y, z, t), s 23 = s 32 = 
P(x, y, z,t), B 23 = B 32 , -B 22 = B 33 , U 23 = U 32 and -IT 22 = n 33 . It 
follows from this and the fact that (T( 2 ) = cr (3 ) that (gTTBJ), (|ODD and flOS]) 
(and hence (|4.19| ) and ( |4.27|) ) are identically satisfied. 

Calculation of (|4.39|) shows that a, (3 must satisfy the Cauchy-Riemann 
equations, 

(4.40) 

ay a z 



d (3 da 
dy d z 







We can write these equations economically as 







(4.4i; 



(4.42) 



with C = y + iz. In addition, combining fl£I5D , (^20|) , (g2| ) and fl07|) 
we find that Il afc and B ab must also satisfy the Cauchy-Riemann equations. 
Thus we have, 

<9n 33 dn 23 

, (4.43) 



dy 



dz 
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dy 


+ 


d z 


dB 33 




dB 23 


dy 




d z 


OB 23 




dB 33 




+ 





du 23 on 33 , , 

(AAA) 

uy u z 

and, 

a R 33 F) R 23 

, (4.45) 

. (4.46) 
uy u z 

This appearance of the Cauchy-Riemann equations is to be expected when 
working with shear-free null hypersurfaces |§, 

Next we examine the integrability condition ( |4.9| ). To evaluate this equa- 
tion we first need to calculate w V c s ab and s ab +6 s ab . The only non-vanishing 
Christoffel symbols for the line element ( |2.1| ) are r{ 4 = A/A, T 2 4 = B/B , 
T 3 4 = C/C, = A A, V\ 2 = BB and Tl 3 = CC. Using these and the 
fact that B = C we find that 

(3) V^ afe = 0, (3) V 4 s afc = 0, (4.47) 

where (3 takes values 1, 2, 3 and 

Sab + 0s ab = ^ + ^s ab . (4.48) 
at A 

Then Eq. ( |4.9|) yields the following: 

with T(t) introduced in (2.20), and 

+S23 (6 1 J 3 + 6 1 b 6 3 a )} , (4.50) 



A 



A 



Hub 



1 ds. 



ab 



A dz 
+s 33 {5151 + 5151)}. 



Since A cab is orthogonal to u c on all its indices we also have, 
Aiab = , A c4b = , A ca4: = . 



(4.5i; 



(4.52) 



With these equations for A cab it is easily shown that Eqs. ( |4.11| ), ( |4. 1 7| ) and 
the last of ( ^.37| ) are identically satisfied. 
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We now turn our attention to ( 4.22 ) and (|4.23 ). We first note that, 



s\ {3) V b A m = 



and 



( 3 ) V a s ab = . 

Then ( |4.23|) reduces to, 

(-^m q Ag m) ^b A q ^mb & m Sqb) 

For convenience we shall define, 



(b)mq 



A 6 — A ° 

■f^m q -^q m 



-X, 



(b)qm ■ 



This allows us to write (|4.55|) as, 

X(b)mq ^b + ( ® q S-mb @ m ^qb) 



(4.53) 
(4.54) 

(4.55) 
(4.56) 
(4.57) 



After a simple calculation using ( 4.49| )-( 4.52j) we find that the non-zero com- 
ponents of X( b ) mq are X^ 12 = 



X(3)i3 and X( 2 )i3 



- ^(3)12 with 



X, 



(2)12 



B 2 



-AB 22 + j A j (A s 22 ) - A (a {2) - a (1) ) s 22 



(4.58) 



and 



X, 



(2)13 



-Ll^-AB 23 + j (J^ + J^) (As 23 )-A(a (2) -a (1) )s 2 ^ , 

(4.59) 

As a result of this X^ mq A& = and ( 4.57|) reduces to 



0" g "Smfc m ^qb . 



(4.60) 



This is clearly an identity since a (2) = C(3) and s 22 = — S33, S23 = S32 are the 
only non-zero components of s a &. Thus we have now shown that (|4.23| ) is 
identically satisfied. That ( |4.22j ) is also identically satisfied can be seen as 
follows: B ab and U ab satisfy the same equations as s ab and so we have, 



a TT _ ^CL TT ri 

" m '-'-aq " q ^^am " 



r-d r> _ a r> n 

V m t^aq ° q *->am u • 



Then Eq. ( |4.22j ) becomes 



(3) V b (A 



— A " \ 

m q -"-q m) 



(4.61) 
(4.62) 

(4.63) 
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which we write as 



X( c ) mq - C + T mc X( c ) qd + T c dc X( d ) mq + Y qc X( c ) dm — , (4.64) 

where X( c ) mq is given by ( |4.56| ). Substituting for X( c ) mq from ( |4.58| )- fl4.59[ ) 
and using the Christoffel symbols which were listed after ( [4.46 ) we find that 
this equation is an identity. 

At this point the equations we have left to satisfy are ( f4.4| ), ( |4.29| ), and 
(|4.42j )- (|4.46| ). As we mentioned previously (|4.29| ) is extremely complicated 
and to simplify the calculations which follow it is convenient to make use of 
the null tetrad in the background space-time which was introduced in ( |2.22|) . 
In terms of this tetrad s ab can be written (because s a a = 0, s ab u b = 0, 
s ab k b = and so s ab l b = 0) 

s ab = sm a m b + sm a m b , (4.65) 

with 

s = -B 2 (a + i(3) . (4.66) 
Since Il a& and B ab also satisfy the conditions given immediately before (|4.65|) 

n ab = Tlm a m b + Tlm a m b , (4.67) 



we can write, 
with 
and 
where 



n = -fi 2 (n 33 + in 23 ) , (4.68) 

B ab = B m a m b + B m a m b , (4.69) 



B = — B 2 (B 33 + i B 23 ) . (4.70) 

We recall that a ab m a m b , which was first introduced in ( |4.34|) , is a gauge- 
invariant first order small variable and we can write it as 

a ab m a m b = KF(<j ) ) , (4.71) 



for some function k. Multiplying the second of (|4. 1| ) by m a m b and replacing 
a ab m a m b by the right-hand side of ( |4.71|) we find (after noting that fn a = 
since m a ]d = B B^ 1 u a m d and m a is orthogonal to w a )that 

k = s , (4.72) 

and 

B = s + 9s. (4.73) 
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We now have an equation for B in terms of s but for clarity we shall refrain 
from replacing B in the equations which follow by the right-hand side of this 
equation until the end. Repeating the above calculation on the first of Q4.1p 
yields 

A cab m a m b = ^V c s , (4.74) 

which is identically satisfied as a consequence of ( |4.49|) -( [4.52D and (|4.73|) . 

We also need to express E ab and H ab on this tetrad. It follows from (|4.21|) 
and the expressions for s ab , B ab and Tl ab given above that 

E ab = Em a m b + Erh a m b , (4.75) 

where 

E= (-U-B) F+ (\s) F' . (4.76) 



.A 

We note in passing that, making use of ( |4. 71| )— ( ^7T2|) and (|4.74| )-( f4T75|) , we 
can write E' ab in ( |4.36|) as 

E' ab = E' m a m b + E' m a m b , (4. 77) 

with 

E'=(~fi-B + js^F+jsF'. (4.78) 

Now we consider H ab . First, substituting for A ca b from Eqs. (f4.49|) -( ^.52|) in 
( |4.12| ) we observe that H ab = except for 

*° = - fa - *») •» + *" " ^ + I) « ->} F 

-]-s 23 F' } (4.79) 

and 

*> = = fa - *») •» " " ^ + ^) F 

+j S 22 F', (4.80) 

where d/d T = A(d/dt). Then in terms of the null tetrad we can write 

H ab = Hm a m b + H m a m b , (4.81) 
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with 

-i-sF' . (4.82) 

We note that since s ab , B ab and IT afe satisfy the Cauchy-Riemann equations, 
s, B and IT are independent of (. But H and are derived from these 
quantities and thus they are also independent of (. 

Next we wish to express W bt given by ( |3.28| ) on the null tetrad but first 
we need to examine h d b i]^ rsd u r H a S]d . Using (|4.81j) and recalling that m a ]d = 
B B~ l u a m d we deduce that 

h b a rf rsd u T H a s . d = r] trsd u r H 4 m b m s + rf rsd u r H 4 m b m s . (4.83) 

Taking into account that dH/d( = we find that, with respect to the null 
tetrad, we can write 

H, d = V2 B- 1 H- C rh d - A' 1 (H x + H T ) u d + H x <j) !d , (4.84) 

where = dH/d(, H x = d H/dx and Ht = d H/dT. Putting this into 
fl£53D yields 

hi V trsd u r H\, d = j (H x m b m l - H x m b m 4 ) , (4.85) 

which is symmetric in b and t. We are now in a position to write W bt in 
terms of the null tetrad. We find it convenient to introduce a new variable 
S which is defined by 

E bt + -U bt = £fh b fh t + £m b m t . (4.86) 

It is easily seen (using ( |4.67| ) and ( |4.76| ) ) that 

S = {U-B)F + A- 1 sF' . (4.87) 

Then replacing E bt + \li bt and h d b r]^ rsd u r H a s , d in (|3~28|) by (|185D and (|486D 
respectively gives (since ( 3 ^V c m a = and m a = 0) 

W bt = (-i--9S-iA- 1 H x ^m b m t + £m f m {t a b)f 

+ (-£ ~^0S + IA- 1 H x ^j m b m l + £ rh f m {t a b)f . (4.88) 
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But a simple calculation shows that 



and thus we can write 



07 2 ) m b m* 



£ _ 1 
B 3 



9 \ m b m t 



where 



= Wm b m t + Wm b m t 



AB dt ' A 



Substituting from (|4.82|) and ( |4.87|) reveals, 

W = PF + QF' 



(4.89) 

(4.90) 
(4.91) 

(4.92) 



with 



1 c) 1 

— -{AB (-B + fi)} + - ((7(2) - (7 (1) ) s a 



+ A Bx A 3 B 2 \<)r or 



(4.93) 



and 



i 



,0-(2) - 0"(1)J S 



A 3 £ 2 \<9x <9T 
1 9 



(4.94) 



We now turn our attention to Eq. (f4.29|) . We first observe that since 
°"(2) = °"(3)) C(i) + 2o"(2) = 0, and s 22 = — s 33 , s 32 = s 23 are the only non- 
zero components of s ab , the coefficient of F' in Eq. (|4.29 ) vanishes identically. 
Therefore 



WV c W bt = 3(-a 2 A c bt + a b f a f s A c st + a b f a st Aj & 
W s cJ st A b f -h bt a a d a af A cfd )F. 



(4.95) 



Using fl4.49|) -( |i.52j) in this we find, after a lengthy calculation, that the co- 
efficient of F on the right hand side vanishes and so we are left with 



(3) v c iy M = o . 

Thus from ( |4.90|) (since ^V c m a = 0) we have 

^w c w = , 



(4.96) 



(4.97) 
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which implies that 



dW 
dx 







dW 







dW 







(4.98) 



Substitution of W from (l492p into the first of (|498|) and using the arbitrari- 
ness of -F(0) we find that we must have P = and Q = 0. Hence W = 
and the other two equations in (|4.98|) are identically satisfied. We find it 
convenient now to write, 



(4.99) 



where Q is an analytic function of (. This follows from ( |4.42j ) and ( [4.661 ). 
Then from ( |4.94| ) with Q = we obtain an equation for S, namely, 



fl = -^( D g + A§ 6 



(4.100) 



Here the operator D is given by D = d/dx + Ad/dt = d/dx + d/dT. It 
follows from this and ([4 . 6 8|) that n 33 + i U 23 is analytic in ( and so ([143]) 
and ([4.44j) are automatically satisfied. 

We now consider ( |4.93|) with P = 0. Substituting for fl and B from ([4.73|) 
and ( |4. 1 00|) respectively results in the following wave-equation forQ(( ,x ,t): 



D 2 Q + A 



B 
B 



A' 
A 



DQ = . 



(4.101) 



We note that we have made use of the field equations (|2 . 3|) — ( |2~5|) (with B = 
C) to write the equation in this form. 

The only equation remaining to be satisfied is ( [4.4] ). We shall first examine 
the case when a = 1. Using ( ^4.49|) - (f4.52| ) and recalling the only non-zero 
Christoffel symbols for the line-element (|2.1| ) with B = C given prior to 
Eq. Q4.47D above, we find that ( ^4.4|) with a = 1 reduces to 



0(2) - o- {1) ] 







(4.102) 
BB- X m a u d 



Multiplying this equation by m a m b and noting that m a ^ 
yields 

(a {2) -a {1) )Dg = . (4.103) 

It is easily checked that Eq. ( |4.103|) is equivalent to Eq. ( f4.102[) . We are in- 
terested here in an anisotropic universe and for this reason we cannot have 
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C(2) = C(i) (we already have a^) = o"(3) and if we also put <7( 2 ) = o"(i) the 
result is an isotropic universe). Thus ( |4.103| ) gives us a remarkably simple 
first order wave-equation, 

DG = 0, (4.104) 

from which we can conclude that Q = Q(( ,x — T(t)). Now ( |4.101| ) is au- 
tomatically satisfied. We note that if the background cosmology is isotropic 
then the line-element (2.1) is simplified to A = B = C = t 2 ^ 3 = R(t) (say). 



In this case a( 2 ) = o"(i) and so (|4.103|) is satisfied while the wave-equation for 
Q(( ,x ,t) is ( |4.101| ) specialised to 



D 2 Q = , (4.105) 

and ( |4.1UU| ) becomes 



ti = -^{Dg + kg) . (4.106) 

These are precisely the equations (5.37) and (5.36) found in [jl]] when our gen- 
eral formulation for studying shear-free gravitational radiation in isotropic 
cosmologies is specialised to apply to perturbations of the Einstein-de Sitter 
universe. 

We have yet to examine ( f4.4| ) when a = 2, 3 or 4. The case when a = 4 



is identically satisfied as a consequence of (^4.52|) . To check that ( [14|) is 
satisfied when a = 2 or a = 3 we first multiply across by m a m b . Noting that 
( 3 )V c m a = this yields 

(3) V C £ = A cab m a m b -2a f b A fac m a m b + 2a f b A afc m a m b 

+a f c A fab m a m b + \e A cab m a m b . (4.107) 

with B given by Eq.( |4.70] ). Then using 

j-(A cab m a m b ) = (^A m - m b -2^A cab m a m b , (4.108) 



we find after a straightforward calculation that fl4.107| ) is identically satisfied. 

Now that s ab , Tl ab and B ab are known we can calculate E' and H and 
form the gauge-invariant electric and magnetic parts of the perturbed Weyl 
tensor which are given by ( |4.71| ) and (|4.81| ) respectively. We can write the 
result compactly as 

EU + iff* = -2 { (A _ JL\ g + _L A. ( g F) \ ma mi , (4.109) 
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The scalar product of this complex-valued tensor with itself clearly van- 
ishes and this is the algebraic property one expects the gauge-invariant part 
of the perturbed Weyl tensor to have if the perturbations describe gravita- 
tional waves [cf. Eq.(5.42) in QXJ to which Eq.( |4.109| ) specialises when the 
background cosmology is Einstein-de Sitter]. 



5 Discussion 

The feature that makes this study of shear-free gravitational radiation 
propagating through a Bianchi I anisotropic universe so much more com- 
plicated than the corresponding study when the universe is isotropic [|T| is 
the unavailability to us of the full perturbed shear of the matter world-lines 
as an Ellis-Bruni variable. This is in turn due to the central role played 
by the perturbed shear of the matter world-lines along with the perturbed 
anisotropic stress of the matter distribution in the study of gravitational 
waves. Nevertheless this paper demonstrates that these difficulties can be 
circumvented and the final result, which is summarised in the perturbed 
anisotropic stress ( [4.1 00| ) and the gauge-invariant perturbation of the Weyl 



tensor ( |4.107| ), with the function Q involved in both satisfying ( (4.1011 ) and 
( |4.102| ), is easily surveyable. These perturbed quantities have been calcu- 
lated with the aid of the basic Ellis-Bruni variables listed in ( |4.1|) . This is 
because A ca b , G ca b , B a b , C a b and Ii a b appearing in ( |4.1| ) are now all known 
in terms of Q: with s afe given by ( |4.65|) and ( f4.99|) we have C a b and G ca b given 
by Q and (pp . In addition B ab is given by fl4J3|) and flOSD while 



A ca b is given by (CT)-(ra) and Ii ab is obtained from (^BTQ and flOUOl) . 
At this level the most significant difference between the perturbations of the 
anisotropic universe compared to those of the isotropic universe is that in 
the former case the final wave equation to be satisfied by Q is first order 
whereas it is second order in the latter case. It is clear from ( (4.103[ ) that 
the anisotropy in the shear of the matter world-lines of the background is 
directly responsible for this. 



Finally we see from ( |2.23| ) that the null hypersurfaces = constant in 



the anisotropic background have null geodesic generators with shear provided 
(T(2) 7^ o"(3) • If is informative to carry out an analysis in parallel with that of 
section IV in this case. The calculations are briefly outlined in Appendix B. 
The significant conclusion is that the gauge-invariant part of the perturbed 
Weyl tensor does not possess all of the algebraic properties associated with 
gravitational waves having propagation direction <p^ a in the anisotropic back- 
ground. For example, we cannot have H ab 0^ vanishing if the null hypersur- 
faces <f)(x a ) = constant we are working with have shear. Therefore we cannot 
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interpret the perturbations as describing gravitational waves in this case. 
This observation suggests to us the possibility of constructing a non-vacuum 
generalisation of the Goldberg-Kerr theorem |TD| which will also generalise 



the theorem of Szekeres |11] expressing the non-existence of Petrov Type N 



dust-filled universes. This will be discussed on another occasion. 
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A Notation and Useful Equations 



Throughout this paper we use the notation and sign conventions of and we 
take the speed of light to be unity. We consider a four dimensional space-time 
manifold with metric tensor components g ab , in a local coordinate system 
{x a }, and a preferred congruence of world-lines tangent to the unit time- 
like vector field with components u a (with u a u a = — 1). With respect to this 
4-velocity field the energy-momentum-stress tensor T ab can be decomposed 
as 

T ab = fiu a u b + ph ab + q a u b + q b u a + n ab , (A.l) 

where 

h ab = g ab + u a u b , (A.2) 

is the projection tensor and 

q a u a = 0, n ab u a = 0, n a a = 0. (A.3) 

Then fi is interpreted as the total energy density measured by an observer 
with 4-velocity u a , q a is the energy flow (such as heat flow) measured by 
this observer, p is the isotropic pressure and 7T a b is the trace-free anisotropic 
stress (due, for example, to viscosity). After absorbing the coupling constant 
into the energy-momentum-stress tensor Einstein's field equations can be 
written as 

Rah ~ ^ 9abR = T ab , (A.4) 

where R ab = R a c bc are the components of the Ricci tensor and R is the Ricci 
scalar. 

We indicate covariant differentiation with a semicolon, partial differentia- 
tion by a comma and covariant differentiation in the direction of u a by a dot. 
As usual square brackets denote skew-symmetrisation and round brackets 
denote symmetrisation. Thus the 4-acceleration of the time-like congruence 
is 

u a :=u a . b u b . (A.5) 
With respect to u a and h ab , u a . b can be decomposed into 

u a -,b = u ab + a ab + ^ 6 h ab - u a u b , (A.6) 

where 

Wat := U[a;b] + U[ a U b ] , (A.7) 

is the vorticity tensor of the congruence, 

&ab ■= U(a;b) + «(a«t) , (A.8) 
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is the shear tensor of the congruence and 

■= u a , a , (A.9) 

is the expansion (or contraction) of the congruence. 

Finally in this section we give the equations obtained when the Bianchi 
identities, written compactly as 

C abcd d = R c[a;b] _ }_ g c[a R ;b] ^ ^ 

are projected along and orthogonal to u a . They are the div-E equation, 
h b g E° d .j h f d + 3co s H b - rf a ™ u a a d p H qd = ±h b c /r c 
+\ {-* M A + u b a cd rr cd - 3 co bd q d + - bd q d ~6q b + n bd u d } ,(A.ll) 

the div-H equation, 



h b g H° d .j h f d -3u s E b + V ba ™ u a a d p E qd = (fi+p) uo b 
+\r,\acU q q a > c + \r l b qac u«{u d « + G dc )K a d , (A. 12) 

the E -equation, 

h) E f9 ^ + h a b rf> sd u r H a s . d - 2 H ib rf^ drs u d u r 

_ E (t w b)s _ 3 E (t a b)s + h tb E dp ^ + E bt = _K fi+p ^ a tb 

- 1 - h tb {fi + 6(n + p)} - q^ b ««> - \ «<> q t] - \ 
A {uJ 4b + a c(b } u t) ^ + K u (t g b) _i^t + n c(p u t) ^ 

2 6 2 

-\{u< b + a< b }^ c -^9^ bt , (A.13) 

and the H -equation, 

h) H f9 h\ - h {b r]^ rsd u r E a s . d + 2 Ef rf^ drs u d u r 

-Hf u b)s - 3 a b)s + h tb H dp a dp + 9 H bt = -q (t cu b) 

-\ V {t rad {oo b)d + a b)d } u r q a - ^ b rad rr^ d u r 

+ ^V {b ra d u t) u r {uj cd + a cd }n a c . (A.14) 
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B Null Hypersurfaces with Shear 



We note from (|2.23| ) that the null hypersurfaces <fi = constant have null 
geodesic generators with shear provided oy 2 ) °"(3)- We briefly describe in 
this appendix the perturbations analogous to those of Section IV in this case. 
Equations (|4.1|) -( P2|) still apply and we shall again assume that B ab tfi^ = 
(the reason for this is given before ( |4.37| )). Now recalling that s ab is orthog- 
onal to u a and trace-free it follows from (|4.10|) and (|4.18| ) that s ab has only 
one independent component s 23 = s 32 . B ah and IT afc are also orthogonal to 
u a and trace-free and satisfy the same equations as s ab . Thus they also have 
only one independent component, B 23 = B 32 and II 23 = n 32 . On account of 
this, the Cauchy-Riemann equations ( |4.40| )-( [4.46| ) now read, 



ds 



23 



dy 
dU 23 

dy 
8B 23 

dy 







ds 



2:1 



dz 
dU 23 

dz 
dB 23 

dz 











(B.l) 



:b.2) 



;b.s) 



In the present context (ct( 2 ) 7^ <7(3)) 
derived from ([4.9|) are modified to: 



.4 



la6 
A2ab 



(KK + KK){-ab 



2.3 



equations ( [4.49|) -( ^:.52|) which we 



^23 , j A ds 23 

+ A s 2 3 + A 



dx 



dt 



As 23 {a {3) -a {1) ){6 l a 6 s + 8l5 l b ) 
As 23 (a {2) -a (1) )(5 1 a 5 2 b +5 2 a 5 1 b ) 



(B.4) 

(B.5) 
(B.6) 

Aiab = , A a4b = , A aM = . (B.7) 

It follows immediately from these equations that ( j4.11|) and (}4.17|) are iden- 
tically satisfied. Using these equations in ( }4.12| )-( [4.14| ) we find that H a b = 
except for 



and 



H 



22 




:b.s) 



H 



3.3 



'AC 
AB 



Ti s 23 



— B03 

B 



B 
A^C 



C 
A 2 B 



D{As 23 )\F 



B 

AC 



C 



(B.9) 



D{As 23 )\F + — s 23 F' 



AB 



(B.10) 
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Here the operator D is given by D = d/dx + Ad/dt. We note that as a 
consequence of these equations we have 

H ab <j>, b = 5l S -f(a {3) -a {2) )F^0. (B.ll) 

Thus one of the algebraic properties of the gauge-invariant part of the 
perturbed Weyl tensor, namely if afe 0,ft = , which we expect to hold if 
the perturbations describe gravitational waves having the null hypersurfaces 
= constant as the histories of their wave-fronts, is not possible unless the 
null hypersurfaces are shear-free. 

We now consider ( [3.281 ). First, a straightforward but tedious calculation 
gives, 

h a(b Vt) rsd u r H\, d = \ A ( A # 33 _ °_ h 22 ^ {6 l 5 l + si 5 2 ) . (B.12) 

Then substituting for E ab and H a t, from ( |4.21|) and (B8)-(B10) respectively 
in ( |3.28| ) and making use of ( |B.12| ) yields 

W ht = (PF + QF')(6 2 b 5* + 6 3 b 5 2 ) , (B.13) 

where (after recalling that + ct( 2 ) + o"(3) = 0) 

d /TT „ x f 1 (B C\ A) /TT „ , ldB 23 

p = - s (n- - fl.) + 1- + _j - - j ( n 23 - b 23 ) + --^ 
-^^^aTj-SI^ to"' (R14) 

and 

1 1 ds 23 1 (B C\ 1 , N 

Q = A^-A^ + 2A\B + c) S ^-A^ D{As ^ 

-2A a V 8m -Ai-te- (R15) 
Remembering that the only non-zero components of the background cr a b 
are = o x x, o"(2) = o" 2 2 and = cr 3 3 , and substituting for A cab from 
(B4)-(B7) we find that, in this case, the right-hand side of ( [4.29 ) vanishes 
identically. Thus we have 

(3) V c H^ = 0. (B.16) 

It is easily shown, using the Cauchy-Riemann equations and the non-zero 
Christoffel symbols given before ( |4.47| ) that this equation is identically sat- 
isfied when c = 2, 3 or 4. However when c = 1 it implies, 

^-o. ,B,r) 
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It follows from this and ( |B.13| ) that P = and Q — 0. Similar to the 
shear-free case we now find it convenient to define 



S23 



and 



B 



23 



Then from ( |B.15|) with Q = we obtain 

n 23 = 



2 nn IB c 



(B.18) 
(B.19) 

(B.20) 



Finally we examine Eq. (O) when the null geodesic generators of the 



hypersurfaces have shear (i.e. when am 7^ C(3))- When c = 2 or c = 3 we 



find, after a lengthy calculation utilising (B4)-(B7), that 7i in ( B.19| ) is now 
given by 



H = 2G t 



A 



Q, 




:B.2r 



where Q t = dQ/dt and Q x = dQ/dx. We note that this equation is valid if 
and only if <j(i) 7^ a (2) and 7^ (T( 3 ). Then substituting for y4 m6 and 7i 
from (B4)-(B7) and (p.21| ) respectively in (jO]) with c = 1 yields 



A 2 



d 2 q 




c 



og 



'AB AC 
AB + AC 








(B.22) 



where we have used the field equations 



of Q. Eqs. (|B.20|) and (|B.21|) allow us to write, 



n 23 — B 



23 



to simplify the coefficient 



^B.23) 



and 



]_dB^__^ ,J_(B C\ 

A dx " A 2 A3^ + A 2 lB + ci^ 



;b.24) 



Using these two equations it is a simple calculation to show that (B14) with 
P = is identically satisfied. We note that ( [4.4|) with c = 4 gives = 0. 

The basic Ellis-Bruni variables listed in ( |4.1| ) are now determined in terms 
of Q, with Q(x, t) given by ( B.22|) .The only non-vanishing component of s a b 
is given by ( |B.18| ) and of B a b by ( |B.19| ) and ( |B.21| ), with the corresponding 
component of U a b given by ( |B.20|) . As in section IV, C a b and G ca b are obtained 
from ( f4.5| ) and ( f4.6| ) while A ca b is given by ( [B.4} )-( |B~7D in this case. 
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